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ABSTRACT 

We present an one-dimensional numerical study of Alfven waves propagating along a 
radial magnetic field. Neglecting losses, any spherical Alfven wave, no matter how small 
its initial amplitude is, becomes nonlinear at sufficiently large radii. From previous 
simulations of Alfven waves in plane parallel atmospheres we did expect the waves 
to steepen and produce current sheets in the nonlinear regime, which was confirmed 
by our new calculations. On the other hand we did find that even the least nonlinear 
waves were damped out almost completely before 10 Rq. A damping of that kind 
is required by models of Alfven wave-driven winds from old low-mass stars as these 
winds are mainly accelerated within a few stellar radii. 
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1 INTRODUCTION 

Spherically symmetric models are often the simplest realistic 
models conceivable for astrophysical systems, and have been 
applied successfully to describe not only hydrostatic objects 
such as stars, but also hydrodynamic processes, for instance 
stellar winds (Parker 1958). In many cases it is desirable to 
extend a spherically symmetric hydrodynamic problem to 
the analogous magnetohydrodynamic problem. In a global 
sense such a magnetohydrodynamic analogue cannot exist, 
as that would require the existence of magnetic monopoles. 
Locally, however, the magnetic field is divergence-free, and 
presents a natural approximation of a region with a diverg- 
ing magnetic field. In hydrodynamics there are only acoustic 
waves (for spherical sound waves see e.g. Landau & Lifshitz 
1987), but the magnetic field introduces Alfven waves. These 
transverse magnetic field oscillations cannot be extended to 
cover an entire spherical surface without introducing discon- 
tinuities, but once again it will not present any difficulties 
locally. 

In a previous paper (Boynton & Torkelsson 1996, here- 
after Paper 1) we have shown that nonlinear Alfven waves 
in a planar geometry can steepen and form current sheets, 
and thereby be damped by Joule dissipation or by doing 
mechanical work on the background medium. The efficiency 
of this mechanism is limited by the Alfven wave becoming 
less nonlinear as it propagates upwards through a strati- 
fied medium. A spherical Alfven wave on the other hand 
may become less nonlinear for some time, but eventually 
has to grow nonlinear again because of the divergence of 
the background magnetic field. Furthermore the gas pres- 
sure decreases faster than the magnetic pressure with height 
independently of whether the symmetry is plane-parallel or 



spherical, so that the Alfven wave becomes more dynami- 
cally important. Models of Alfven wave-driven outfiows have 
been constructed by Hartmann & MacGregor (1980, 1982) 
to explain the winds from late-type giants, although a more 
popular model is that the outflows are driven by radiation 
pressure working on dust in the stellar atmosphere. It has 
been pointed out that the Alfven waves must damp within a 
few stellar radii to avoid accelerating the wind to too high ve- 
locities. There is no generally accepted model for this damp- 
ing and Holzer, Fla & Leer (1983) have pointed out that the 
stellar wind depends sensitively on the damping mechanism. 

A spherically symmetric model may also apply to the 
coronal holes observed on the Sun (Bohlin 1976, Zirker 
1977). The magnetic field in a coronal hole is open, which al- 
lows the plasma to expand outwards and form the so-called 
fast solar wind. This component of the solar wind is too fast 
to be described by the cleissical Parker (1958) model, and it 
is possible that the extra acceleration is provided by Alfven 
waves. Alfven waves have been observed further out in the 
solar wind (Belcher & Davis 1971, Balogh et al. 1995), but 
it is not known how these waves relate to the ones that are 
supposedly present in the corona. 

Section 2 describes the initial hydrostatic model and 
summarizes the properties of linear waves propagating 
through it. The results of our numerical simulations are pre- 
sented in Sect. 3. Section 4 discusses the physical interpre- 
tation and the astrophysical consequences of our results, in 
particular, with respect to solar and stellar winds, and our 
conclusions are summarized in Sect. 5. 
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2 THE STATIC MODEL AND ITS WAVE 
MODES 

2.1 The static model 

We assume a spherically symmetric Sun-like star. On top of 
the stellar surface there is an isothermal corona at a tem- 
perature of 10^ K. For a corona in hydrostatic equilibrium 
the density can be written as 

R z 



po (z) = po (0) exp ^- 



H R + z 



(1) 



where po (0) is the coronal density at the stellar surface, R 
is the stellar radius, z the height above the stellar surface 
and H = 2kBTR^ /{GMmn) = 6.110'' m, a scale height, 
with kB Boltzmann's constant, T the temperature, G the 
gravitational constant, M the mass of the star, and mn 
the mass of a hydrogen atom. Throughout this paper sub- 
script Os, with a few obvious exceptions, refer to the undis- 
turbed background and Os as arguments to functions mean 
the functional values at z = 0, that is the stellar surface. 
The isothermal sound speed, Cs, which is what matters as 
we are assuming isothermality in the dynamical model, is 
1.3 lO'^ ms"^ and po(0) is put to 510"" kgm"^. A ver- 
tical magnetic field, Bz{z) oc [R + z)~^ , gives an Alfven 
speed, va = Bz/^/iop, with the z-dependence of Fig. 
For -Bz(O) — 3 10~^ T, the Alfven wave is supersonic even 
at z = URq, whereas for 5^(0) = 10"'^ T the Alfven 
wave is subsonic everywhere. Assuming a period of 300 s 
the wavelength is always small compared to the scale height 
(Fig. §b), except for 5^(0) = 3 10"" T when it exceeds the 
scale height for z < 2Rq and is comparable to the solar 
radius at z « 1.5i?Q. Note that for 5^(0) = 10"^ T we 
choose a period of 900 s, to get the same wavelength as for 
Bz(0) = 310"^ T. 



2.2 Linear theory of spherical waves 

The magnetohydrodynamic equations for a spherically sym- 
metric system can be written as 

dp 1 d 



dt (R + zf dz 



[{R + zf pv,] =0, 



d 



1 d 

{R + zf dz 



1 d 



^'^'-pitrz' 



[(7? -I- z)^ P'VyVz'] = 

')) -71 _ 

-B^J^ - p- 



R + z' 



dp 



--B.Jy + ByJ.-pg + p-^., 



1 d 



R + zdz 



[{R + z) B,v, 



(2) 



(3) 



(4) 



(5) 



(6) 
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Figure 1. (a) Alfven velocity as a function of height, z, for mag- 
netic field strengths of 310—* T (solid line), 3 10-^ T (dashed 
line) and 1 10"^ T (dotted line) as measured at z = 0. The sound 
speed is plotted as a dash-dotted line, (b) The wavelength of an 
Alfven wave with a period of 300 s as a function of z for the 
same magnetic field strengths. We do not plot the wavelength for 
1 10-5 T, as its period of 900 s gives it the same wavelength as 
310-5 T with a period of 300 s. The dash-dotted line denotes 
the local pressure scale height and the long-dashed line the solar 
radius 



dBy 



+ 



d 



dt R + zdz' 
1 d 



R + zdz 



{R + Z) ByVz] = 

[{R+z)B.Vy]., 



(7) 



where p, Vx, Vy, Bx and By stand for the density and the 
transverse components of the velocity and the magnetic 
field. Vz is the velocity along the background magnetic field, 
g the gravitational field strength, and p = c^p is the pressure. 
The need for an energy equation is eliminated by assuming 
all processes to be isothermal, and B^ must be independent 
of time to keep the magnetic field divergence-free. The elec- 
tric current density is given by 



tJ X 



and 



d 



po R + z dz 



[{R + Z) By 



d 



po R + z dz 



[(R + z) Bx], 



(8) 



(9) 



where po is the permeability of free space. 

We linearize the equations around a, possibly, stratified 
static background medium and a radial magnetic field B^ — 
Bz{0) ("h^) • For brevity we leave out the y-components, 
and use r — R + z. 
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Figure 2. The maximal period for acoustic waves in a stratified 
atmosphere, Ps, as a function of z 



dp I d , 2 \ „ 



dt 
d 



at \r J HQ r dr 



dt 



2 dp 

-Cs-^ P9, 

or 



dB^ 1 d f fR\^ \ 



(10) 

(11) 
(12) 

(13) 



where we have written the density as p = po + p. These 
equations separate into two groups, Eqs. (Mand (^ de- 
scribing acoustic waves, and Eqs. (^l|) and (1^) describing 
Alfven waves. 

As a first example we assume that po is constant and 
g = 0. The acoustic waves are then described by two wave 
equations 

d'p 



dt^ 
and 



_l_d_ 

^2 



^ 2d_ 

dt^ dr 



1 d 



J.2 

which have the solutions 

~ ~ /r\\ i(kr — ujt) 

p = p{0)—e^ 



(14) 



(15) 



(16) 



where p(0) is the amplitude of the density fluctuations at 
r — R and 



LJ H — 

po dr 



^i(kr — ujt) 



ik^r 

with the dispersion relation 

2 2,2 



(17) 



(18) 



(cf. Landau & Lifshitz 1987, Ch. 70). 

In the case of a stratified medium the wave equation for 
the sound waves can be written as 

d^p cl d 



d P ^ S^d_ (^2dp\ ^ dp 

Q^2 j,2 Qj- \ Qj- J Qf 



(19) 



For g = GM/r^ the background density is oc exp (— -^-h^); 
where H is the scale height at z = 0. We write the density 
fluctuations as 



10.00 
^~ 1.00 
05 0.10 
0.01 



10 12 14 



1000 


00 ! 


100 


00 r 

\ 


10 


00 r 


1 


00 r 





10 





01 - 




10 12 



14 



10000.0 
1000.0 






c 


100.0 
^ 10.0 








1.0 








0.1 









4 6 8 10 
2 (solar radii) 



12 



14 



Figure 3. (a) The relative amplitude of an Alfven wave prop- 
agating through a stratified atmosphere. The solid, dashed and 
dotted lines represent waves of initial relative amplitudes of 1, 
0.1 and 0.01, respectively, (b, c) The ratio of gas to magnetic 
pressure for the background field, (b) and the oscillating field 
Sj: (c). is 310""' T (solid Une), 3 lO'^ T (dashed line), and 
dotted line. The relative amplitude of the oscillations is 
B,(0)/B,(0) = 0.1 



p = p(0):^exp(-A_^) 



i(fcz-cjt) 



2H R + z, 

(cf. Paper 1). The resulting dispersion relation is 

clR^ 



2,2 , 



4/^2^4 ' 



(20) 



(21) 



implying that acoustic waves with frequencies u < Ns — 
CsR^/{2Hr^) are evanescent (cf. Lamb 1908, 1932). We plot 
the correspondiiig period, Ps = 27v/Ns in Fig. 

Equations (hll) and ( |l3| ) can be rewritten as (cf. Heine- 
marm & Olbert 1980, Leer, Holzer & Fla 1982, MacGregor 
& Charbonneau 1994) 



dl 
dt 

and 

dl ^ 

dt " 



VA 



■ VA 



<»'(f)( 
'»)(7)'( 



dl 
dr 



1 + 
r dr 



dg 



(22) 



(23) 
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Table 1. Simulations of linearly polarized Alfven waves with period, P, travelling through an isothermal medium stratified according 
to Eq. (|^) with a base density 5 10~^^ kgm""^ and a background magnetic field, Bz , decreasing as 1/{R + z)'^ . The temperature of the 
medium is 10® K, which corresponds to a sound speed of 1.2910^ ms~^. The computational domain covers the the corona between 1 
and 15 Rq outside of a solarlike star for Models la - Ic, and between 1 and 16 Rq for Models 2a - 3c. The homogeneous (H) models are 
similar but has a constant density of 10~^^ kgm~^, their radial extents are 0.15 Rq (Hla-Hlc), 0.3 Rq (Hid) and 0.6 Rq (H2a-H2c) 
with the lower boundary at IRq. The driven Alfven wave is polarized at 45° to the x-axes for all of the simulations. At denotes the 
time step and N the number of spatial grid points used in the respective models. (0) is the strength of the background magnetic field, 
VA (0) the Alfven velocity, (-Bosc / Bz ) (0) the relative amplitude of the Alfven wave and /3mag (0) the plasma beta all calculated at 2 = 



Model At (s) N B,{0) (T) P (s) VAiO) (ms-i) {Bosc/B.){0) Anag(O) 



Hla 


0.25 


720 


510-6 


300 


4.410* 


0.01 


17 


Hlb 


0.25 


720 


510-6 


300 


4.410* 


0.1 


17 


Hlc 


0.25 


720 


510-6 


300 


4.410* 


1.0 


17 


Hid 


0.25 


1440 


510-6 


300 


4.410* 


1.0 


17 


H2a 


1 


720 


2 10--"^ 


300 


1.8105 


0.01 


1.0 


H2b 


1 


720 


210-5 


300 


1.8105 


0.1 


1.0 


H2c 


1 


720 


210-5 


300 


1.8105 


1.0 


1.0 


la 


1 


15120 


1 10-5 


900 


1.310* 


0.01 


210 


lb 


1 


15120 


1 10-5 


900 


1.310* 


0.1 


210 


Ic 


1 


15120 


1 10-5 


900 


1.310* 


1.0 


210 


2a 


1 


15120 


310-5 


300 


3.810* 


0.01 


23 


2b 


1 


15120 


310-5 


300 


3.810* 


0.1 


23 


2c 


1 


15120 


310-5 


300 


3.810* 


1.0 


23 


3a 


0.25 


7560 


310-"' 


300 


3.8105 


0.01 


0.23 


3b 


0.25 


7560 


310-"^ 


300 


3.8105 


0.1 


0.23 


3c 


0.25 


7560 


310-"' 


300 


3.8105 


1.0 


0.23 



where 

/ = Wo; 



VMoPo 



(24) 



and va(0) — Bz (0) /\//iopo (0) is the Alfven velocity at r = 
R, so that 



- if + g) 

and 



(25) 



(26) 



Bx = 2 V^'opo i~f + g) ■ 

Note that f — and g = yield waves propagating 
downwards and upwards, respectively. Substituting /(r, t) = 



F(r)e-'"' and g{r,t) = G(r)e- 



we get 



_d_ 

dr 



F 



G 



F 



G 



(27) 



The general solution to these equations is rather compli- 
cated, and we will restrict ourselves to cite some general 
and useful results, that can be derived from simple physical 
arguments. 

First we derive the scaling properties of the amplitudes 
of Vx and Bx in the WKB sense. Consider the Poynting flux 

S^^, (28) 
Mo 

with the electric field, E = —v x B + J /a, where we assume 
that the conductivity, a, is infinite. If there are no losses, 
the amplitude of the radial component of the Poynting vec- 
tor Sz will be proportional to Sz is also proportional 



to BzBxVx, the oscillating transverse magnetic and velocity 
fields, and as Bz oc l/r'^ in our monopole geometry, BxVx 
is not explicitly dependent upon r. The oscillating magnetic 
and velocity fields of an outgoing Alfven wave are related by 



Bx 



(29) 



which gives Bx oc pj^* and Vx oc Pq We plot the relative 
amplitude of the Alfven wave, Bx/Bz, in Fig. |a. In Figs. |b 
and c we compare Bz and Bx, respectively, to the gas pres- 
sure by calculating the magnetic beta /3x,z = 2^op{z) / Bx^z- 
The main effect of the deviations from the WKB- 
approximation is that the waves can be refiected against 
gradients in the Alfven speed. This problem has been dis- 
cussed for linearized spherical Alfven waves by An et al. 
(1990) and Lou & Rosner (1994). The refiection is important 
as it increases the momentum transferred from the Alfven 
wave to the medium. 



3 RESULTS OF THE NUMERICAL 
SIMULATIONS 

We modified the code of Paper 1 to simulate Alfven waves 
in a spherically symmetric magnetic field. In addition we 
now solve for both transverse components of the velocity 
and magnetic fields, so that also circularly polarized waves 
can be described. In general our waves are linearly polar- 
ized at 45° to the x-axis. Initial test runs showed that the 
boundary conditions from Paper 1 were too reflective at the 
upper boundary. After some experiments we found that a 
satisfactory solution was to extrapolate the boundary val- 
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Figure 4. Magnetohydrodynamic waves propagating in a spheri- 
cally symmetric magnetic field with a uniform background density 
(Model Hid), (a) Ap (solid line), and (dotted line) oscillations 
at 6 000 s. Both quantities are measured in velocity units. The 
lower two panels show Bx (solid line) and Vx (dotted line) oscilla- 
tions at (b) 6 000 s and (c) 10 000 s, respectively, both measured 
in velocity units, (a) shows the acoustic precursor {z > 0.1/Jq) 
to the Alfven wave in (b), and the oscillations generated by the 
magnetic pressure (z < O.li?0). The Alfven wave is gradually 
damped as it steepens into current sheets (c) 



ues of pvz, Bx and By from the last three grid points, and 
calculate pv^ and pVy from the requirement that they should 
describe Alfven waves propagating out through the bound- 
ary. The boundary conditions at the lower boundary are the 
same as in Paper 1, that is we drive an upwards propagat- 
ing Alfven wave, and keep all other variables fixed at their 
initial values. 

Table |^ describes the models that we have calculated. 
The H-models are homogeneous in the sense that they have 
constant density and temperature, and no gravity, but the 
magnetic field still goes as l/(i? -f . 
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Figure 5. Magnetohydrodynamic waves propagating in a spheri- 
cally symmetric magnetic field with a uniform background density 
(Model H2b). The left axis and the solid line show Ap/pg and the 
right axis and the dotted line Bx- The times are (a) 400 s, (b) 
1200 s, (c) 1800 s, and (d) 3 000 s, respectively. At the start 
the Alfven wave is propagating faster than a sound wave so that 
in (a) we see density oscillations carried by the magnetic field, 
but as the Alfven velocity decreases there is a sound wave that 
overtakes the Alfven wave (b-d) 



Table 2. The maximum of the precursory density fluctuations 
(Ap/po)iiiiax ^^'^ second-order fluctuations {^p/po)2min 

in the interacting Alfven and acoustic waves for Models Hla-Hlc 



(Bosc/Bz)(0) 


(Ap/po)lmax 


(Ap/P0)2min 


0.01 


4.810-"^ 


-610-^ 


0.1 


4.8 10""' 


-610-"* 


1.0 


4.710-2 


-610-2 



3.1 Wave propagation 
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Both the precursors and the density fluctuations inside the 
Alfven wave are proportional to the square of the amplitude 
of the Alfven wave. This is natural as they are produced by 
the magnetic pressure of the Alfven wave. 

The situation is more complicated in Models H2a-c, 
where «a > Cs at z = 0, but Cs > va for z > 0.04i?Q. 
The acoustic wave lags behind the Alfven wave at first. The 
density oscillations propagating with the Alfven wave in Fig 
^ are not acoustic waves, but rather p and oscillations 
caused by the magnetic pressure oscillations (cf. HoUweg 
1971). Density maxima occur initially where there is a max- 
imum in the magnetic pressure, but drift out of phase when 
the acoustic wave overtakes the Alfven wave (Fig. |b,c). The 
density fluctuations simultaneously increase in amplitude, 
and eventually steepen into shocks and lose energy (Fig. 
^). It is instructive to study these oscillations as a function 
of the amplitude of the Alfven waves. For a low amplitude 
Alfven wave the density fiuctuations are essentially a propa- 
gating wave which does not transport any significant amount 
of mass (Fig. ^). This wave is growing in amplitude up to 
z = Q.3Rq after which it remains constant or even decreases 
in amplitude again (Fig. ^b). The amplitude of the wave is 
varying roughly as the square of the amplitude of the Alfven 
wave as can be found by comparing Figs, ^a and b, and it is 
propagating at the sound speed, 1.310^ ms~^. pv^ changes 
character completely for a nonlinear Alfven wave (Fig. pc). 
There is an outfiow of mass, which is concentrated to shocks 
that are separated in time by roughly half the period of the 
Alfven wave. The front of the outfiow is initially propagating 
at 2 lO'^ ms"\ but slows down to 1.5 lO'' ms"\ at 3 000 s, 
which is still supersonic. However the outfiow velocity, Vz , is 
mainly subsonic, although it may reach peak values as high 
as 1.410^ ms-\ 



Figure 6. pv^ as a function of z for Models H2a - c (a - c) at 

2 400 s. Note how pvz changes character from being a propagating 
wave (a,b) to an outflow (c) when the amplitude of the Alfven 
wave increases 

3.1.1 Unstratified models 

The amplitudes of the transverse velocities and magnetic 
fields are not affected by the divergence of the radial mag- 
netic field, but they are affected by the stratification in the 
same way as in the plane-parallel case considered in Paper 
1. The most straightforward cases are obviously the ones 
without stratification, where the Alfven velocity decreases as 
l/{R + z)^. Models Hla-d are particularly simple as va < Cs 
everywhere in these models, so that the Alfven wave has an 
acoustic precursor. For high amplitude waves, such as Hlc 
and Hid, both the acoustic precursor and the Alfven wave 
steepen to form shocks and current sheets, respectively (Fig. 

Note that we have transformed the density, Ap = p — po, 
and magnetic field, Bx, oscillations to velocities by multi- 
plying with Cs/po and 1/^poPo, respectively. We quantify 
the amplitude of the precursor by defining ( — ) as the 

y Po ' Imax 

maximum of the density fiuctuations in the acoustic precur- 
sor (the density does not drop below its background value) . 
Likewise we call the minimum of the density fluctuations 
inside the Alfven wave (■^) . ■ The numerical values for 
these quantities are given in I for Models Hla - Hlc. 



3.1.2 Stratified models 

The stratified models are different with respect to the Alfven 
velocity in the sense that va first increases with z, but starts 
to decrease again above z = '^Rq ■ Models la-c have va < Cs 
everywhere, so that we expect to see an acoustic precur- 
sor. We increased the period to 900 s for these models, as 
a shorter wavelength would have required a finer resolution, 
and thus more CPU time. We show snapshots of Models lb 
and c in Fig. ^. Figure ^ shows clear signs of the nonlinear 
damping beyond z — IRq as Vx stops increasing in am- 
plitude. The front of the Alfven wave has advanced further 
in Fig. ^ than ^b, which is surprising as both snapshots 
are taken after the same simulation time. It is the wave in 
Fig. ^ that is propagating at the expected Alfven speed, 
whereas t he n onlinear Alfven wave is propagating too fast 
(see Sect. 

Models 2a-2c are more interesting because the Alfven 
speed is comparable to the sound speed in the interval 
1^0 < 2 < SRq. The most striking feature in the early 
stages of Model 2b is however unrelated to this. Figure ^ 
shows how the acoustic oscillations described by pv^ grow 
in amplitude, but reach a maximum at z = 0.7. After this 
the sound waves steepen and dissipate the momentum that 
they are carrying (Fig. ^). At 50 000 s we see that the ampli- 
tude of the Alfven wave decreases sharply between z = 3Rq 
and z = 5i?Q (Fig. |b). There are also signs of oscillations 
at higher frequencies superposed on the wave, in particular 
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Figure 7. Vz and at 50 000 s for Models lb (a) and (b), and Ic (c) and (d). (a, b) show an Alfven wave of low amplitude with an 
acoustic precursor. At higher amplitude (c, d) the Alfven wave is strongly damped by nonlinear steepening, and the Alfven wave front 
is propagating faster 



downward of z = 5Rq, and a general low- frequency mod- 
ulation at z < IRq. Fig. ^ suggests that the first 7 or 
so wavelengths of the wave are less damped and show less 
extraneous oscillations than the following wavelengths. The 
same effects are present in Model 2c too, but in addition 
the wavefront of Model 2c is propagating super- Alfvenically, 
which we will try to explain in Sect. 4.2, Model 2c crashed 



eventually as the magnetic pressure evacuated a part of the 
grid, which was expected as the pressure of the oscillatory 
magnetic field is comparable to the gas pressure. 

The only one of Models 3 that has got {B^ + By)/{2iio) 
significantly weaker than the gas pressure is Model 3a. There 
are no acoustic precursors to the Alfven waves in these mod- 
els as va > Cs everywhere, and indeed the most advanced 
parts of the density oscillations are not sound waves, but 
density fluctuations carried by the Alfven wave at the Alfven 
speed. Such a density fiuctuation is expected to obey the re- 
lationship 



Ap 

P 



Vz 
VA 



(30) 



(HoUweg 1971, Eq. (15)). At 2 500 s this relationship is sat- 
isfied for z > 2.5-Rq, but at smaller heights there is a large 
density excess (Fig. [To[ ) . The density enhancement grows in 
extent with time (Fig. |l^). The local maximum is moving 
outwards at the sound speed, 1.3 10^ ms^^. All the Models 
3 crash eventually. Model 3a and b at the time that they hit 
the upper boundary, and Model 3c at a much earlier time, 
when the Alfven wave evacuates a part of the grid. 



3.2 Energetics 

As in Paper 1 the equations for the magnetic and kinetic 
energies are 



and 

§i " -^-^P+^-iJ ^ B)+pv-g,{32) 

where S = E x B/no is the Poynting vector, and the elec- 
tric field is given by E = —v x B + J /a. The Ohmic dissipa- 
tion J^/cr in Eq. represents the effect of the numerical 
diffusion in our code. The numerical diffusion smears out 
discontinuities over the length scale of the grid spacing, but 
it is at the same time negligible in regions with smooth ve- 
locities and magnetic fields. The thickness of a current sheet 
is proportional to I/ct, from which it follows that J oc cr, so 
that the dissipation integrated over the current sheet is in- 
dependent of a. The magnitude of the numerical diffusion 
should thus not affect the net damping unless the diffusion is 
strong enough to damp out smooth variations, which is not 
the case. As was pointed out in Paper 1 it is impossible to 
calculate the Ohmic dissipation directly from the quantities 
of the numerical simulation, however the other terms in Eq. 
( ^l[ ) are accessible. We average these terms over the period 
of the Alfven wave. The averages of the nonlinear models 
are sensitive to secular changes and aperiodic fluctuations 
in the simulations, which introduce some uncertainties in 
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Figure 8. pv^ (solid line/left scale) and Bx (dotted line/right 
scale) for the Alfven wave of Model 2b. (a), (b) and (c) are for 
the times 7200, 8 800 and 11200 s, respectively, (d) Shows an 
enlargement of the region between O.5f?0 and IRq in (c). Note 
how the acoustic oscillation (solid line) increases in amplitude up 
until z = 0.7i?Q, after which it decreases again, possibly due to 
dissipation in shocks (d) 



the calculated energy losses. We plot the time average of 
Model 2b as an illustrative example in Fig. 10 - 20 % 
of the Poynting flux is dissipated immediately at the lower 
boundary due to imperfections in the boundary conditions 
and initial state, and it may therefore be advisable to think 
of the waves as being of correspondingly lower amplitude 
than indicated from Tab. 1. Most of the remaining Poynt- 
ing flux is lost between 3 and 5 Rq by doing work on the 
background medium via the Lorentz force. The damping be- 
comes more efficient with increasing amplitude of the Alfven 
wave both in the sense that more of the flux is lost from the 
wave and in the sense that the damping sets in earlier (Tab. 
b|). In addition a larger part of the Poynting flux is spent on 
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Figure 9. Bx as a function of z for Model 2b at 50 000 s. (a) 
< z < 3Rq, (b) 3 < z < 6Rq and (c) 6 < z < 9Rq. The 
Alfven wave steepens to a square wave (a), and loses energy in 
the current sheets (b), but the first few wavelengths in the head 
are essentially unaffected (c). Below z = IRq, in particular, there 
are signs of a low-frequency modulation of the Alfven wave 



doing mechanical work instead of being lost through Joule 
dissipation as the amplitude increases. 



4 DISCUSSION 

The most important and interesting results from our simu- 
lations are 

• The fact that practically all our Alfven waves lose a 
signiflcant fraction of their Poynting flux within less than 
10 Rq 

• The fact that a nonlinear Alfven wave can propagate 
super- Alfvenically 

• The appearance of oscillations at lower frequencies than 
that at which the Alfven wave is driven 

We will discuss each of these facts below, and also the sig- 
nificance of our results for stellar wind models. 



4.1 The physical mechanism of the wave damping 

It was suggested in Paper 1 that the damping of the Alfven 
waves takes place in current sheets. These current sheets ap- 
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Table 3. A compilation of the energetics of the Alfven waves. For every Model the Poynting flux is given as the equivalent Poynting flux 
at the stellar surface, Sz ( ) , S'z(O) is the Poynting flux at 2: = 0, and (final) the Poynting flux after the damping of the Alfven 
wave. We also give the work done by the wave via the Lorentz force IVlof per unit surface area and time, and the position where half of 
the damping has taken place, Zdamp 



Model 


(Wm-2) 


S^final) {Wm-2) 


VV-Lor (Wm-2) 




la 


410-^5 


310-5 


0.410-5 


8.2 


lb 


410-3 


310-3 


110-3 


7.2 


Ic 


410-1 


1 10-1 


310-1 


1.4 


2a 


910-* 


0.410-'' 


410-4 


9.3 


2b 


910-2 


0.510-2 


410-2 


4.4 


2c 


9 


.2 


9 


0.3 





3x10"* 


'In 


2x10^ 




1x10* 









-Ixio"^ 


H 


-2x10* 




-3x10* 





























'1 






A 




/\ 




1 \ 










\' > 












'1 






w 




( 1 




1 1 










'\ / 


1 

I / 








1 


1 


1 1 
















1 1 


' aJ 




1 2 3 4 5 6 

2 (solar radii) 

Figure 10. A snapshot of (a) Bx (solid line) and Vx (dashed 
line), and (b) Vz (solid line) and Ap (dashed line) at 2 500 s for 
Model 3a. Note that that Bx and Ap are measured in velocity 
units by multiplying with 1/y'popO and Cs/po, respectively, (a) 
shows the typical behaviour of an Alfven wave, while (b) shows 
density fluctuations carried by the Alfven wave at ^ > 2.5i?0 and 
an outflow at smaller z 
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Figure 11. Ap/po at 2 500 (soUd Une), 7 500 (dashed Une) and 
12 500 s (dot-dashed line) in Model 3a. The peak of Ap/p is mov- 
ing at the isothermal sound speed 




pear at the nodes of the Alfven wave, because the nodes rep- 
resent minima of magnetic pressure, and thus oppositely di- 
rected magnetic field lines are pushed together at the nodes. 
This mechanism is still at work in our new simulations, but 
at the same time, due to the stratification of the medium, 
there is a gradient in the wave pressure, which in itself may 
do work on the medium (cf. Jacques 1977). In the nonlinear 
regime these effects are coupled together and cannot be sep- 
arated easily. Table ^ shows that the Alfven wave damping is 
a nonlinear effect as the fraction of the Poynting ffux that is 
lost increases with increasing initial Poynting flux. In partic- 
ular it is interesting to find that Whoi increases faster than 
Sz{0), which we interpret as the transverse magnetic field 



Figure 12. Sz ( ^j;" ) (solid line) and the integral of the work 
done by the Lorentz force (dashed line) averaged over 300 s (the 
period of the Alfven wave) for Model 2b 

becoming more dynamically important compared to the gas 
pressure. The fact that Zdamp decreases with ^^(O) shows 
that the damping cannot be described by a linear model. 

4.2 Super-AIfvenic motion 

We have seen in some of our simulations that nonlinear 
Alfven waves can in some circumstances propagate faster 
than the local Alfven speed (Fig. ^a-c). The reason that 
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Figure 13. for Models 2a-c at 50 000 s (a - c), and for 
Model 2c (d). Note the change of Vx-sca\e between (b) and (c), 
and the change of z-scale between (c) and (d). This Fig. shows 
that the wave front propagates faster the higher the amplitude of 
the wave 



the wave in Model 2c can propagate at a higher speed than 
those in Models 2a and b is that the wave itself is acceler- 
ating an outflow with a maximum velocity of 4.8 10'* ms~^ 
(Fig. Thus the Alfven wave gains the velocity of the 

medium it is propagating in. This effect was not found in 
Paper 1, as in the plane-parallel models the Alfven velocity 
is increasing upwards to eventually become supersonic. Con- 



sequently the Alfven wave front manages to stay in front of 
the main outflow (cf. Paper 1, Fig. 7), which is generated by 
the nonhnear part of the wave at lower altitudes. 



4.3 Low-frequency modulations 

As already noticed there is some evidence in Fig. ^ that the 
Alfven wave is modulated on a frequency lower than the 
driving frequency. To test this we calculate power spectra of 
Models lb and 2b (Fig. 0). In addition to the main peak at 
0.003 2 Hz we see side peaks at 0.003 and 0.000 6 Hz. To 
explain these extra peaks we note that a wave reflected back 
down by an inhomogeneity propagating upwards is Doppler- 
shifted to a lower frequency 



VA 



VA + Vz 



(33) 



where v (0) is the original frequency of the wave. To explain 
the side peak in Model lb we require an outflow velocity 

= 0.04wa, which agrees well with the velocities gener- 
ated in the model (Fig. ^). The peak at 0.000 6 Hz of 
Model 2b on the other hand requires — 0.7wAiWhich is 
significantly higher than what is available (Fig. Il4d). An 
alternative interpretation in this case is that the frequency 
0.000 6 Hz represents the beat frequency between the orig- 
inal and reflected Alfven waves, thus giving the frequency 
of the reflected Alfven wave as 0.002 6 Hz, corresponding to 

= O.lvA- This agrees well with the velocity before the 
jump in Hz at 11 Rq (Fig. [l4|d), but the peak in the power 
spectrum (Fig. is at 0.002 Hz, which may be a conse- 
quence of the inhomogeneity of the medium. 



4.4 Importance for stellar winds 

Alfven-wave driven winds have been discussed several times 
in the past (e.g. Hartmann & MacGregor 1980, 1982; Leer 
et al. 1982, MacGregor & Charbonneau 1994). The main 
reason is that the Alfven waves may be able to drive the 
wind in circumstances where more well-understood mecha- 
nisms can be shown to be insufficient. This is for instance 
the case for the fast solar wind moving at velocities of 800 
kms~^, which is too fast for the classical thermally driven 
wind model by Parker (1958), and the winds of late- type gi- 
ants, where the temperatures are too low to drive a wind. A 
severe weakness of the models of the winds of late- type stars 
is that the Alfven waves must be damped within a few stellar 
radii in order to avoid over-accelerating the wind. Usually 
this has been done by ascribing an arbitrary damping length 
to the Alfven waves, but the resulting models are sensitive 
to the choice of the damping length (Holzer et al. 1983). In 
our simulations the waves are damped without the action 
of any dissipative mechanism while they are doing work on 
the background medium. Our mechanism does not require 
the wavelength to become comparable to the length scale of 
the Alfven velocity variations in contrast to some previous 
attempts that have appealed to non-WKB waves being re- 
flected. Although we have not surveyed the parameter space 
it appears likely that the parameters of our model can be 
modified to provide a sufficiently effective damping to ful- 
fill the constraints put on those wind models. An obvious 
improvement on our current simulations is to start from an 
initial state which describes an outflowing wind. 
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Figure 14. Power spectra of Bx for Models lb (a) and 2b (c) at t = 100 000 s. (b) and (d) show in units of the local Alfven velocity 
at the same time for Models lb and 2b, respectively 



Magnetohydrodynamic waves may have been observed 
in the inner parts of the solar wind. Ofman & Davila (1997b) 
have calculated that magnetohydrodynamic waves in the so- 
lar wind may give a line broadening of ~ 300 kms"'^, which 
is comparable to the line widths observed by SOHO (e.g. 
Kohl et al. 1996). In particular Ofman & Davila (1997a) 
suggests that what they call solitary waves are essential for 
the acceleration of the solar wind. The solitary waves are 
propagating p-Wz -perturbations generated by Alfven waves, 
and thus are akin to the density oscillations we have found in 
for instance Fig. ^. The authors call them solitary waves as 
they believe them to be related to the solitary wave solutions 
found in slabs and thin flux tubes by Roberts (1981) and 
Roberts & Mangeney (1982). This is an interpretation that 
we find questionable for the following reason. It is certainly 
true that the two-dimensional axisymmetric model used by 
Ofman & Davila bears some resemblance to a flux tube with 
a lower density in the interior, but our one-dimensional sim- 
ulations lack all features characteristic of a flux tube, and 
still we flnd the same kind of density oscillations propagat- 
ing faster than the sound speed. We thus conclude that the 
oscillations cannot be related to the existence of solitons in 
flux tubes. Regarding the question of whether the solar wind 
is driven by Alfven waves or solitary waves that is very much 
a matter of semantics as the solitary waves must be driven 
by the Alfven waves. It is interesting to compare the sim- 
ulations in terms of efficiency of deriving energy from the 



Alfven waves. The background models are roughly similar 
apart from that Ofman & Davila assume a radial magnetic 
field which is 70% stronger than the one we use in Models 
3. They drive the Alfven wave at an amplitude inbetween 
our Models 3a and b, and assume a period which is almost 
an order of magnitude longer than ours. Consequently they 
find larger radial velocities than we do, but unfortunately 
they do not give any numbers for the fraction of the Poynt- 
ing flux that has been converted to mechanical energy. They 
do however state that 3 10~^ of the Poynting flux has been 
lost due to Ohmic dissipation. This appears to be a rather 
inefficient conversion process, but one must keep in mind 
that their grid extends only to 4 Rq, and we typically find 
that most of the damping takes place at larger distances. 

Naturally there are alternative models for driving the 
fast solar wind. Feldman et al. (1996) have suggested that 
the outflows are driven by the same kind of reconnection 
events that produce the X-ray jets (Yokoyama & Shibata 
1995). In this model the emerging magnetic field of a bipolar 
region collides with the magnetic field of the chromospheric 
network. This field is concentrated in narrow fiux tubes in 
the photosphere, but due to the stratification the flux tube 
increases in radius upwards. At the collision point the mag- 
netic fields reconnect and two jets moving at close to the 
Alfven velocity appear. Only a part of the momentum of 
the jet is directed upwards, but the horizontal momentum is 
absorbed by nearby matter and magnetic field, so that also 
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a part of this becomes available for generating an outflow 
from the coronal hole. Inevitably some of the energy goes 
into producing Alfven waves, and so the model may also 
explain the presence of Alfven waves in the solax wind. 



5 CONCLUSIONS 

In this paper we have studied the propagation of nonlinear 

spherical Alfven waves. Like in our previous simulations of 
Alfven waves in a plane-parallel atmosphere the waves damp 
by forming current sheets in which Poynting flux is lost to 
Ohmic heating and the acceleration of an outflow. In general 
most of the Poynting flux is spent on accelerating an out- 
flow. This combined process of damping Alfven waves and 
accelerating an outflow may be important in understanding 
both the fast solar wind and the winds of late-type giants. 
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